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Abstract 

Acoustic attenuation measurements in air have been analyzed in order to estimate the second 
coefficient of viscosity. Data over a temperature range of 11°C to 50°C and at relative humidities 
between 6 percent and 91 percent were used. This analysis has shown that the second coefficient 
of viscosity varied between 1900 and 20,000 times larger than the dynamic or first coefficient of 
viscosity over the temperature and humidity range of the data. In addition, the data have shown 
that the molecular relaxation effects, which are responsible for the magnitude of the second 
coefficient of viscosity, place severe limits on the use of time-independent, thermodynamic 
equations of state. Compressible flows containing large streamwise velocity gradients, like shock 
waves, which cause significant changes in particle properties to occur during time intervals shorter 
than hundredths of seconds. must be modeled using dynamic equations of state. The dynamic 
model approach has been described briefly. 
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Introduction 

When the stress tensor is symmetric, continuum theory limits the number of arbitrary 
constitutive constants to two for a linear, isotropic material. While solid mechanics handbooks 
tabulate pairs of material constants (Young’s modulus and Poisson’s ratio) relating stress to strain, 
fluid mechanics handbooks tabulate only one constant — dynamic viscosity - - relating stress to 
rate of strain. Because many important fluid flows can either be treated as incompressible flows 
or boundary layer flows, the influence of a second coefficient of viscosity on the equations of 
motion can be eliminated from the formulation and the problem is avoided. However, current 
research for a variety of vehicle systems in air require analyses of flow fields which are neither 
incompressible nor amenable to boundary layer approximations. While limited experimental data 
are available, the purpose of this work is to examine recent acoustical data in order to estimate 
the magnitude of the second coefficient of viscosity in air. 

Since the second coefficient of viscosity appears as a linear combination with the dynamic 
viscosity in many tensor manipulations, the linear combination is often referred to as the bulk 
viscosity or volume viscosity . Furthermore, since the trace of the stress tensor is an invariant, 
an average normal stress or pressure force is predicted for most compressible fluid flows which 
occurs in addition to thermodynamic pressure. That problem is avoided when the ratio of the 
second coefficient of viscosity to the dynamic viscosity is negative two thirds. In 1845, C.G. 
Stokes suggested that using a second coefficient equal to -2/3/i was a convenient way to avoid 
uncertainties in the meaning of pressure in a flowing fluid. Stokes' hypothesis can be justified 
for monatomic molecules like helium, but no such arguments exist for mixtures of polyatomic 
molecules like air. 

White 1 has provided a good review of the second coefficient of viscosity problem in fluid 
mechanics. Since normal shock waves produce large velocity gradients in the direction of flow, 
shock thickness measurements can be used to determine the influence of second coefficient of 
viscosity. Unfortunately, shock waves produce strong gradients in all of the thermodynamic 
variables and if the second coefficient of viscosity varies with temperature, the accuracy of those 
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estimates are limited. However, the early shock wave thickness measurements of Sherman 
which employed thermocouples and position indicators for the measurements, showed that the 
second coefficient of viscosity in air was a positive value (not equal to -2/3/i). 

Bulk viscous effects in air are of major concern to acousticians. Attenuation of radiated 
sound waves is controlled for the most part by bulk viscosity. However, because frequency 
dependent thermophysical properties violate the fundamental laws of continuum mechanics, the 
acoustics literature is nearly void of references to bulk viscous effects. Acousticians use total 
sound attenuation (in nepers/meter, where nepers are dimensionless) and tabulate that property. 
Total sound attenuation in air varies strongly with acoustic frequency, relative humidity and 
temperature. As productive engineers, acousticians have been concerned more with propagation 
models than with the frequency dependent property inconsistencies. 

The second coefficient of viscosity controversy may have reached maximum intensity in 
1954, when Rosenhead 3 convened a meeting of some of the leading authorities in continuum 
theory and fluid mechanics in an attempt to sort out the problem. A prominent contributor to 
the controversy has been Truesdell, 4,5 who showed rigorously that bulk viscous effects in simple 
fluids cannot vary with frequency. The stature of the scholars involved in that controversy may 
have restricted the rate at which others contributed information to the subject. Goldstein, 6 was 
willing to discuss the dilemma in 1972, but his publication was rather obscure. 

Acousticians made major strides in unraveling the controversy in the mid 1960’s. Bauer 
and Roesler 7 reported on the slow rate at which mixtures of nitrogen and oxygen returned to 
thermal equilibrium. The reason for slow relaxation rates could be attributed to the vibrational 
degrees of freedom possessed by both oxygen and nitrogen molecules. They found that up 
to 20 million collisions were required to re-equilibrate vibrationally excited molecules. Those 
measurements showed that equilibrium based equations of state could not be used in mixtures 
of nitrogen and oxygen for processes with characteristic frequencies in excess of a few hundred 
herz. Subsequently, Bauer 3 and Kneser* developed rigorous models for these systems using 
ideas from non-equilibrium thermodynamics (see deGroot and Mazur 10 ). While rigorous, the 
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models are cumbersome and demand additional thermodynamic relaxation data that were not 
available for air. 

The implications of this previous work for compressible, non-boundary layer flows in air are 
profound. First, if bulk viscosity fails to obey Stokes hypothesis, a great deal of computational 
analysis is flawed. Second, if molecular relaxation effects are important at frequencies higher 
than a few hundred herz, simple continuum theory models are restricted severely in both space 
(mean-free-path constraints) and time. Finally, attenuation data suggest that the second coefficient 
of viscosity in air may be free of frequency dependent, non-equilibrium relaxation effects only 
at low frequencies. 

Fortunately, the situation in air may not be as bad as it sounds, which is probably why 
the controversy persists. First, as will be discussed later, water vapor accelerates the return 
to equilibrium of nitrogen and oxygen molecules. Hence, below altitudes of 10.5km (35,000 
ft.) water vapor levels are often high enough to increase the collision-driven, upper frequency 
limit to several thousand herz (see Bass et al. 11 ). Second, most computational fluid dynamics 
codes have used artificial viscosity to control computational stability and in some cases those 
adjustments may be consistent with air data. Also, as mentioned previously, incompressible 
flows and boundary layer flows are not influenced by bulk viscosity, although constraints on 
boundary layer approximation limits are more severe. 

Recently, Marcy 12 has reported estimates of the bulk viscosity in air. He employed the 
'correlations of Bass et al. 11 , who used data from Zuckerwar and Meredith 13,14 , to develop 
estimates of bulk viscosity over a wide range of temperatures and relative humidities. This 
study has employed Zuckerwar’s experimental data directly to establish the accuracy of the 
empirically based predictions. 

The present work will start with the basic formulations from acoustics and fluid mechanics, 
for the purpose of unifying the constitutive constants. Using the acoustic and fluid mechanical 
interpretations, Zuckerwar’s measurements have been translated into bulk viscosity data for air. 
Finally, a brief discussion of non-equilibrium thermodynamic effects has been included. 


Formulation 


The stress tensor, <7jj is assumed related to the rate of strain tensor by: 


<7i, = p 


5ui 5uj 5u k 

^ + + 7, W' j 


( 1 .) 


where p is the dynamic viscosity, the u; are the Cartesian velocity components, *-□ is the 
dimensionless form of the second coefficient of viscosity and is the Kronecker delta. The 
Einstein summation convention is used throughout. 


Since the trace of a second rank tensor is invariant, an average pressure contribution, p\ 
due to viscous effects is given by: 


P 


= — crkk/3 = —p{2/3 + 


( 2 .) 


or from conservation of mass, the divergence of velocity is given by: 


<9u k 1 D/j 

9x k p Dt 

so that viscous pressure fluctuations can be represented by 


P ; = v (2/3 + 7n) 


D p 

Dt’ 


(4.) 


where u is the kinematic viscosity and material time derivative of density in an 

Eulerian system. 

The bulk viscosity, f3, is defined by 


0 s p (2/3 + Til) 


(5.) 


and the dimensional form of the second coefficient of viscosity is often designated as A (=/<■• u )• 
Acoustic disturbances are normally very small compared to other fluid mechanical effects. It 
consideration is given to a weak, plane acoustic wave, travelling in the x-direction, the position ot 
an air “particle” translates across a reference location by some amount £(x, t), where the particle 
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translations can be influenced by distance from the acoustic source - - say x. Furthermore, if a 
steady wave train is considered, the conservation of mass can be approximated by 


Po dx 


( 6 .) 


where p 0 is the undisturbed density and p is the instantaneous density. 

If the speed of sound is c, acoustic pressure fluctuations, p, are related to density fluctuations 

by: 


2 a 2 df 

p ss p 0 c A % -p 0 c‘ — . 


(7.) 


The linearized conservation of momentum equation for a viscous fluid obeying Eq. (1.) can 
be written in terms of local particle displacement as 


d 2 £ _ l dp da n _ 2^i, (o . ^ d 2 ( dj 

dt* ~ p 0 dx + dx c dx 2 + + ^ l) dx 2 vat 

where variations of i/ and 7n with temperature have been neglected. 

For an acoustic plane-wave being generated continuously at x=0, 

£ = £ o e- 0,x e i < wt - kx > 


( 8 .) 


(9.) 


satisfies the governing equations where a is called the total sound attenuation (in nepers/meter) 
by acousticians. Substitution of equation (9.) into equations (7.) and (8.) yields 

1 e 


a 

u> c 


' /5 v /(i + * 2 ) + (i + « 2 ) 3/J 


where e = i/(2 + 7n) w /c 2 


( 10 .) 


(H.) 


In most acoustics applications, it is presumed that e is small so that to a good approximation, 
Eq. (10.) becomes 

a i/(2 + 7n) 




2c 3 


( 12 .) 
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When the approximation is met and the viscosities are constant, it is seen that acoustic attenuation 
increases quadratically with frequency. In fact, a great deal of experimental attenuation data are 
plotted as a/f 2 where f is frequency. It should be noted that Stokes derived this relation for 
- ;l l=-2/3 in 1845. 

Air is composed primarily of oxygen and nitrogen with variable amounts of water vapor. 
Both oxygen and nitrogen are diatomic molecules with one vibrational, two rotational and three 
translational degrees of freedom. On the other hand water vapor is triatomic and is classified 
as an angular triatomic molecule 9 . Bass et al. 11 have tabulated the various types of molecular 
relaxation processes which can occur when significant populations of vibrationally excited (out 
of equilibrium) molecules are present in air. Both vibrationally excited oxygen and nitrogen 
molecules relax slowly in air. Surprisingly, water vapor molecules can interact with either 
excited nitrogen or excited oxygen molecules to accelerate their relaxation. Unfortunately, water 
vapor levels of less than 0.01% are insufficient to “catalyze” the relaxation process effectively and 
hence, using the U.S. Standard Atmosphere 15 and calculating maximum water content, dry air 
must be assumed at altitudes above 10,500 m (35,000 ft.). Acousticians have developed empirical 
data to estimate attenuation over a range of temperature, pressure and humidity conditions and 
a typical correlation is discussed subsequently. 

If <j> is used as a measure of the mole fraction of water in air, then 

d > = (Relative Humidity in percent) x P w /Patm (13.) 

where P w is the saturation pressure of water corresponding to the atmospheric temperature 
and Patm is the atmospheric pressure. The influence of water vapor, ambient pressure and 
temperature on total sound attenuation have been correlated empirically by Bass et al. 11 using 
o as an independent variable. When pressures are normalized using atmospheres and absolute 
temperature, in degrees Kelvin, is normalized with respect to 293. 15K, i.e. 


0 = T/'293.15 K, 


(14.) 
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Bass et al. 11 reported total sound attenuation data correlations for air, given by: 


a 


1.83x10 


-11 © 


1/2 


(2x) 


+ © 


-5/2 


0.01278- 


,— 7.638/S 


fr, O + — "jj. 

(2t) fr,o 


+ 0.1069 


,-11.43/e 


fr,x + 




( 2 5r ) 2 fr,N 


(15.) 


where 


fr, o 


P 


24 + 44, 1000 


0.05+0 
0.391 + 0 


(16.) 


and 

fr,N = P (9 + 2000) (17.) 

which are the relaxation frequencies for oxygen and nitrogen, respectively. Hence, if a/w 2 or 
a/f 2 is a constant in an acoustical total sound attenuation plot and e, from Eq. (11), is small, 


71 1 



-2 


(18.) 


Otherwise, Eq. 


(10.) and (15.) can be used to find t and 
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7H = 2 

ijJV 


(19.) 


Equation (15.) has been used to calculate 7n, which has been compared to the experimental 
data. At high altitudes, dry air can be assumed and Marcy 12 has used the data of Bass et al." 
to calculate the second coefficient of viscosity for dry air. 


Reduction of Acoustic Data 

A typical total attenuation plot from Zuckerwar and Meredith 13 is shown in Figure 1. If 
e, defined in Eq. (11.), is small, then from Equation (12.), acoustic attenuation, which is con- 
trolled by simple continuum theory, should vary quadratically with frequency. Furthermore, as 
frequency is increased, vibrational non-equilibrium effects become significant and a simple, time- 
independent, continuum theory cannot be used. Hence, a straight line has been drawn through 
the low frequency portion of the data, in Figure 1 showing how the quadratic approximation 
fits. The appearance of noticeable vibrational non-equilibrium effects are indicated with an ar- 
row where the data cease to follow the quadratic line. From those data, second coefficient of 


Absorption Coefficient, db/m 
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viscosity can be determined. The data in Table I are a summary of the results from Zuckerwar’s 
measurements along with a comparison of estimates of 7n using Eq. (15). 



Frequency, hz 


Figure 1 Typical Attenuation Curve 


TABLE I 

Second Coefficient of Viscosity 
for Air, Calculated from Zuckerwar 13 
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Experimental Conditions 
T P 

Relative 

Humidity 

-%xl0 1<3 

t/xlO 5 

fmax 

711 

fmuXlO 4 

Error in 
Eq. 15 

11.5°C 

1.0009 atm 

9.7% 

37.5 

1.456 m 2 /s 

70 hz 

19,900 

11.1 

-16% 

19.6 

1.0104 

7.9 

14.4 

1.495 

170 

7,770 

9.14 

156 

19.8 

1.0036 

65.3 

5.17 

1.531 

410 

2,730 

10.6 

-71 

20.0 

1.0138 

7.1 

17.6 

1.494 

130 

9,520 

9.88 

134 

20.7 

1.0104 

28.3 

9.73 

1.514 

180 

5,210 

7.56 

-34 

20.9 

1.0172 

73.5 

3.95 

1.526 

420 

2,100 

7.16 

-71 

21.1 

1.0104 

20.4 

11.0 

1.514 

180 

5,900 

8.54 

3 

21.7 

1.0261 

35.5 

7.04 

1.504 

400 

3,820 

12.1 

-76 

21.8 

1.0145 

45.8 

5.83 

1.527 

420 

3,120 

10.6 

-54 

21.8 

0.9458 

91.2 

3.95 

1.661 

480 

1,940 

8.20 

-80 

29.8 

0.9744 

24.2 

12.5 

1.662 

200 

6,390 

11.0 

-41 

29.8 

1.0179 

38.3 

7.04 

1.602 

320 

3,730 

9.88 

-57 

30.0 

1.0118 

52.1 

5.17 

1.625 

410 

2,700 

9.30 

-67 

30.2 

1.0091 

46.0 

5.17 

1.626 

490 

2,700 

11.1 

-58 

30.3 

0.9934 

13.7 

12.5 

1.627 

180 

6.540 

9.88 

43 

30.7 

1.0084 

65.5 

3.95 

1.649 

500 

2,040 

8.68 

-73 

39.2 

0.9989 

25.7 

6.59 

1.725 

480 

3,400 

14.1 

-21 

39.4 

1.0125 

47.9 

5.17 

1.734 

500 

2,650 

11.5 

-70 

39.6 

1.007 

9.4 

20.7 

1.695 

130 

10,900 

12.0 

28 

39.6 

1.0132 

16.5 

10.5 

1.693 

230 

5,520 

10.8 

7 

40.0 

1.0132 

32.9 

5.83 

1.719 

500 

3,030 

13.0 

-46 

49.1 

1.0091 

28.8 

5.17 

1.834 

850 

2,620 

19.9 

-37 

49.8 

1.007 

17.4 

10.1 

1.821 

500 

5,180 

22.8 

-18 

50.2 

1.0234 

6.2 

20.7 

1.769 

180 

11,000 

16.8 

97 


It can be seen from Table I that estimated values for 7 n range from 1940 at 21.8°C, 91.2% 
relative humidity up to 19,900 at 11.5°C, 9.7% relative humidity. The data in Table I are plotted 
as 7 n vs temperature in Figure 2. Because of the extreme influence of moisture content on * u 
at a given temperature, the measured value of relative humidity is indicated beside each data 
point. In addition, the upper, dry-air limit line and lower, saturated air limit line, developed 
by Marcy 12 are shown. While the table shows that Eq. (15) is not very accurate in correlating 
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- n calculated from Zuckerwar and Meredith’s 13 experiments, the variation of 711 with humidity, 
which covers more than two orders of magnitude at a given temperature, can produce much 
larger differences than those resulting from use of Eq. (15). It is important to recognize further 
that the experimental data are self consistent, showing decreases in 711 with increases in water 
content, and that the data fall comfortably within the limit curves predicted by Marcy 12 . 



Figure 2 Variation of non-dimensional second coefficient 
of viscosity with temperature and relative humidity 

When the more subtle frequency limitations indicated in Table I, are included, it should 
be recognized that both the simple time independent continuum model and the implied (from 
Eq. (12)) value of bulk viscosity cease to represent the flow physics fully for a variety of 
unsteady flow conditions. When flow events occur in air with characteristics times on the order 
of milliseconds, molecular relaxation effects must be included directly in the flow model. A 
brief discussion of how relaxation effects can be included follows. 


Thermodynamic Considerations 
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The data discussed thus far suggest that non-equilibrium processes occur during acoustic 
excitations above a few hundred hz in ambient air. Acousticians are concerned primarily 
with higher frequencies and must model relaxation processes directly. Here, non-equilibnum 
is assumed to mean a process which occurs so rapidly that individual molecules do not have 
time to re-establish equilibrium energy distributions over the length scales associated with the 
process. Relying on irreversible thermodynamics, Bauer 8 formulated a dynamic equation of 
state which includes relaxation effects. For the purposes of this discussion, suppose we consider 

a molecule like oxygen which can exist in a rotationally excited state - - say Oo and a 

vibrationally excited state - - O". Then, the following “excitation reactions” are possible due 
to collisions: 

M + 0 2 ^ Oj + M , V12 — -V 2 \ ( 2 ° ) 

M + O2 — O2* + M , V 2 z = -vn ( 21 -) 


M + 0 2 — 0” + M , v\z — -~v$i 


( 22 .) 


where M represents any colliding molecule, and the i/y will be used to designate the degree 
to which each of the oxygen reactions proceeds. Obviously, excitation reactions with other 
molecules are also possible (see Bass et al. 11 ). Furthermore, when significant populations of 
excited molecules are present, thermodynamic temperature ceases to be a unique quantity. In 
order to develop an appropriate model for high frequency acoustic waves, shock waves or other 
flows where bulk viscous effects are important, the second law of thermodynamics must be 
employed. That development has been discussed thoroughly by deGroot and Mazur 10 . They 
showed that the rate of entropy production, T, is given by: 


r = 4 (*) - k t jk « [*£ (f ) - H - | r ' a ' <”• 
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where q is the heat flux vector, Jk is the diffusion flux of the Kth species, given by: 

Jk = PK (uk - u). ( 24 -) 

/ 4 K is the thermodynamic or chemical potential of the Kth species (out of N total species), Fk 
is the body force acting on the Kth species, <Tij is given by Eq. (1.), Ri is the reaction rate (in 
mass production rate per unit volume) of the Ith species, out of M possible reactions, and the 
chemical affinity, Aj, is defined by 

N 

Ai = ^2 v KIPK ( 25 -) 

K=1 

where the vu are the degrees of reaction (indicated in Eqs. (20)-(22». Capital letter subscripts 
do not employ the Einstein summation convention. 

If attention is focused on a single reaction - - say Eq. (20.) - - in a multi-component system, 
and viscous effects are ignored, it is possible to develop a dynamic energy equation model which 
incorporates relaxation effects. Under the single reaction, no viscous effects assumptions, Ri = 
r and all other Ri are ignored. Furthermore, 


N 


Eq. (23) simplifies to: 


A = Ai = ^ vki PK 
k=l 


r = 

T 


(26.) 


with 


where 


dCi 
1 dt 


= r, 


(27.) 


(28.) 


(29.) 


Ci = pi/p and p = 

l=i 

Now the iqi ^ the degrees of reaction from species “I” to “1” and are positive if “1” is a 
product and negative if it is a reactant. It is assumed that the are normalized with respect 
to the molecular weight of I (Mi) such that 

q N N 

Jet'll = - 1 , = 1. andj^i/n = 0, 

1=1 1 I=q+ 1 1=1 


( 30 .) 
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where the vu are an-anged so that reactants are indicated for I < q and products for I > q. 


From Equation (28.) 



dC 1 T 1 9 \r 

ST ' 1 = 1,2 N ' 


we can define a new variable, (, by 


and integration yields 


1 dC\ _ d( _ t 

v\\ dt ~ dt p ’ 


C(t) = 


Ci(t)-C|(0) 

I'll 


(31.) 


(32.) 


We see that for any known Q, all other Q(t) are known when the initial distributions of Q(0) 
are known. Also, we note that any null value for a un results in a requirement that Cj(t) = 
Cj(0). Since Cj < 1 for all t, from the definition Eq. (36.), 


and 


l-Ci(O) <c< _Ci(0) 

i'll ~ “ «di 

Ci(0) c i(°) 

I'll “ ~ I'll 


1 < I < q 
q < I < N 


(33.) 


Now from Gibbs equation, 


N 


(34.) 


de = Tds - pdv + ^ — dCi, 

1=1 A 

where e is the specific internal energy, we can incorporate Eq.s (26), (29) and (31) and 
differentiate with respect to time to write: 


Ds 1 De p Dv A D£ 

Dt = fDt + TDT _ TDt 

Hence, from the exact differential properties of thermodynamic variables, 

>-4Sl.=(SL. 


(35.) 


(36.) 
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(37.) 


dhA _ f _fd g' 

^C/ S ,p \^C/v,T V^/p,T 

using Maxwell relations, where h is enthalpy, a is Helmholz free energy and g is Gibbs free 

energy. It is also noted that A=0 at equilibrium. 

From elementary thermodynamics, an equilibrium can be determined for any gas mixture at a 
prescribed state. Let C eq be the equilibrium value for ( corresponding to any pair of independent 
thermodynamic variables. If we assume that the system is perturbed away from equilibrium 
and that it returns to equilibrium via a constant density (specific volume) and entropy process, 
we can write: 


where A; = ( - ( eq - Hence, from Eq. (31) 


r dQ 
p dt 


and deGroot and Mazur 10 (p 34) have shown that 


M 


r = -£/,jAj/T = 
J=1 


= _/A = jY^e\ eq 
- T -tUV s ,v 


for a constant volume (or temperature) process. As a consequence, 

dC / /d 2 e\ eq ‘ 

at" p7\dc- ) c 


/ s,v 


or defining 


_ 


mz 


(38.) 


(39.) 


(40.) 


(41.) 


(42.) 


we have 


A( = AC(0)e~ t/r *’ v (43.) 

This is an expression for the relaxation time required to re-equilibrate a fluid undergoing a single 
excitation reaction which is not too far from equilibrium. 
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Bauer 8 has shown that the relaxation times of the type displayed in Eq. (42) can be 
incorporated in a dynamic equation of state. If it is assumed that z = z(x,y) where z is the 
dependent thermodynamic variable (e, h, s, g, etc.) and x and y are the independent state 
variables (x = v or - p, y = s or T) and subscripts o are used to indicate equilibrium states, 


1 + r xy 


D 

Dt 


(z - Zo) = 


+ 


dz\ 

dx) 

dz 


( dz 


a/Wc, 7 *"' 


D 

Dt 


% 


A,x 


dz\ 


D 


ir ) + "et: r *y 


fly/c,* 


Dt 


(x - Xo) 


(y - yo) 


(44.) 


Further simplifications can be made for acoustic waves, but they do not apply to shock waves 
or vortices. 

Relaxation time constant measurements which could be used in a dynamic equauon of 
state like Eq. (44) are limited. Bauer and Roesler 7 measured relaxation constants for mixtures 
of nitrogen and oxygen at 300 K. Their measurements showed that the internal temperature, 
Tj, relaxed toward its equilibrium temperature, T , 0 (=T), for an isothermal, constant pressure 
process according to: 


Ti (t) - T io = [Ti(0) - Tije c ? 


(45.) 


where C° is the steady-state, constant-pressure specific heat and C£° is the instantaneous, infinite 
frequency limit, specific heat. For a mixture of 21 percent oxygen and 79 percent nitrogen, their 


measurements suggest that the isothermal response constant was r° = 0.032 s. 

Adiabatic relaxation processes are of more interest than isothermal processes because an 
external energy source is not generally available to maintain isothermal conditions. In that case 
(see Kneser 9 ): 


,DTi 


/p Dt 


(46.) 


-l[c» p (T,-T i0 ) — RT(l — 

where the constant pressure specific heats, C£° and C£ correspond to the instantaneous (infinite 
frequency) and steady-state (zero frequency) specific heats, respectively. Furthermore, if attention 


9 . 


is restricted to an isentropic process, we can write 9 : 


, D 
1 + Tps Dt 


(p ~ Po) = 


dp\ 


3pA, + rp 'l3py<,DtJ 


d£\ D_ 


(P “ Po) 


(47.) 
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where 


Tps = C~r°/C 


(48.) 


Acousticians identify the two density derivatives as equilibrium and infinite frequency sonic 
speed quantities, where 


dp 

dp 


£ ) = 1/c l and 


A,s 



(49.) 


and recognizing that p-p Q is the sound pressure, p, along with incorporating Eq. (6) in Eq. 


(47), they write: 


1 + Tps — 


D_ 

Dt 


P °dx 


1 r ps D 

c 2 c 2 Dt 


P 


(50.) 


They consider the linearized form of Eq. (50), along with the inviscid momentum equation {u= 0 
in Eq. (8)) to be a complete set. Hence, they assume relaxation effects are controlling and 
neglect bulk viscous effects in the conservation of momentum equation. 


Conclusions 


The analysis presented here supports the contention that molecular relaxation processes 
occurring in air result in bulk viscous effects which invalidate Stokes 7 hypothesis. Viscous 
dilatational stresses become important in compressible flows and can easily exceed viscous shear 
stresses in flows with significant streamwise velocity gradients. Due to the strong influence of 
moisture content on bulk viscosity, great care is required in comparing theoretical predictions 
with experiment. 

This work has shown that a serious concern to the aeronautics community should be the 
limited range of compressible flow conditions over which time independent equations of state can 
be employed. The influence of dilatational effects on pressure, along with the slow vibrational 
relaxation rate for oxygen and nitrogen molecules in air, must be accounted for in dynamic 
energy and momentum equation formulations. 
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